Damping of condensate collective modes due to equilibration with the non-condensate 
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We consider the damping of condensate collective modes in the collisionless regime at finite temper- 
atures arising from lack of equilibrium between the condensate and the non-condensate atoms, an 
effect that is ignored in the usual discussion of the collisionless region. As a first approximation, we 
ignore the dynamics of the thermal cloud. Our calculations should be applicable to collective modes 
of the condensate which are oscillating out-of-phase with the thermal cloud. We obtain a generalized 
Stringari equation of motion for the condensate at finite temperatures, which includes a damping 
term associated with the fact that the condensate is not in diffusive equilibrium with the static 
thermal cloud. This inter-component collisional damping of the condensate modes is comparable in 
magnitude to the Landau damping considered in the recent literature. 



PACS numbers(s): 03.75.Fi. 05.30.Jp 67.40.Db 
I. INTRODUCTION 

The collective oscillations of a condensate at zero tem- 
perature T = are well described by the solutions of the 
linearized Gross-Pitaevskii (GP) time-dependent equa- 
tion of motion for the condensate wavefunction $(r, t). 
At finite temperatures, the condensate dynamics is mod- 
ified by interactions with the non-condensate atoms in 
the thermal cloud, which has the effect of renormaliz- 
ing and damping the condensate oscillations. Recently 
the coupled dynamics of the condensate and the thermal 
cloud has been the subject of several theoretical stud- 
ies |Q-§|. Such calculations lead to a generalized GP 
equation for <fr(r,t) and some appropriate Boltzmann- 
like kinetic equation describing the dynamics of the non- 
condensate atoms. In the present work, we make use of 
the recent formulation of Zaremba, Nikuni, and Griffin 
(ZNG) || to discuss a new kind of damping of condensate 
oscillations that arises from collisions between the con- 
densate and non-condensate components. In contrast to 
Ref. S, which discussed the collision-dominated hydro- 
dynamic regime, here we discuss the collisionless regime. 
Within the well-known Thomas-Fermi (TF) approxima- 
tion, we derive a generalized Stringari wave equation de- 
scribing the condensate normal modes B that is valid at 
finite T and includes damping due to the fact that the 
condensate is not in equilibrium with the thermal cloud. 
This new source of damping is in addition to the usual 
Landau and Beliaev damping considered in the collision- 
less region at finite T ]7|-|l4[] . 

Our theory can be used to generalize any discussion 
based on the usual GP equation at T = 0. This simplic- 
ity is due to our neglect of any dynamics of the thermal 
cloud. Available studies of collective modes || at finite 
T suggest that, for any given mode symmetry, one mode 
mainly involves motion of the condensate (with a small 
out-of-phase motion of the non-condensate). This mode 
is a natural extension of the T = oscillation of a pure 
condensate and should be described by our theory. The 



other mode, of the same symmetry, mainly involves the 
motion of the thermal cloud (with a small in-phase mo- 
tion of the condensate) and can be viewed as the natural 
extension of the oscillations above the critical tempera- 
ture Tbec |{l5|l6 |. Our present calculations do not ap- 
ply to such "normal fluid" oscillations, which include the 
Kohn mode at the trap frequency. 

II. DERIVATION OF MODEL 

Our starting point is the finite T generalized GP equa- 
tion derived by ZNG (see also Refs. Q] and ||) 
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Ucxt + gn c + 2gn — ihR 
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where the interaction parameter g = Ai:Ti 2 a/m, a is the 
s-wave scattering length, n c (r,f) = |$(r,i)| 2 , and n(r,t) 
is the non-condensate local density. The damping term 
in (0) is given by R(r,t) = T 12 {r,t)/2n c (r,t), with 



dp 



r Ci 2 [/(p,r,t),$(r,t)]. 



(2) 



This involves the collision integral Cn [f, $] describ- 
ing collisions of condensate atoms with the thermal 
atoms, which also enters the approximate semi-classical 
kinetic equation for the single-particle distribution func- 
tion (valid for k^T 3> gn c o and k^T 3> Tiuiq) 
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V p / = C 12 [/,$] + C 22 [/,<&]. (3) 



Here the collision integral C 22 [/] describes binary colli- 
sions between non-condensate atoms. It does not change 
the number of condensate atoms and hence does not ap- 
pear explicitly in (Q). These coupled equations p)-(|^), 
along with equations (23a) and (23b) of Ref. || defin- 
ing the collision integrals Ci 2 and C 22 , were derived in 
the semi-classical approximation. However, they are ex- 
pected to contain all the essential physics in trapped 
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Bose-condensed gases at finite T, in both the collision- 
less and hydrodynamic domains. They assume that the 
atoms in the thermal cloud are well-described by the 
single-particle Hartree-Fock spectrum e p (r, t) — p 2 /2m + 
U(r,t), where U(r,t) = U exb (r)+2g[n c (r,t)+fi(r,t)}. We 
expect this semi-classical description to break down only 
for very low temperatures where the Bogoliubov excita- 
tion spectrum is more appropriate |l7| ]. Our entire dis- 
cussion is within what is called the Popov approximation 
in that we have ignored all effects associated with the 
anomalous pair correlations m(r,t) = (ip(r, t)"0(r, <)). 

The coupled equations (|l|)-(^) have been used to de- 
rive the generalized two-fluid hydrodynamic equations 
in the collision-dominated region described by a local- 
equilibrium Bose distribution pjl^| . They also have been 
recently used to give a detailed analysis of condensate 
growth by quenching the thermal cloud distribution jfjj . 
In these papers, (l)-(3) are solved with both the conden- 
sate and non-condensate being treated dynamically and 
allowed to be out of equilibrium. The key limitation of 
the present paper is that we only consider the dynam- 
ics of the condensate, with the thermal cloud being in 
static equilibrium. This assumption allows a simple the- 
oretical development and should be adequate for out-of- 
phase modes. The key role of the condensate and non- 
condensate being out of diffusive equilibrium was first 
stressed in a series of papers by Gardiner and cowork- 
ers |[|,^0) ■ These were based on a kinetic master equation 
formalism quite different from what we use, and no appli- 
cation was made to the damping of condensate collective 
modes. 

It is important to understand what is meant by the col- 
lisionless regime and to clarify how this terminology re- 
lates to the present work. Above Tbec (where C12 = 0), 
equation (|^) reduces to the Boltzmann equation describ- 
ing a normal gas 0,0. In this case, the collisionless re- 
gion is well defined and corresponds to having LOiT c \ 3> 1, 
where oji is the collective mode frequency of the gas, on 
the order of the trap frequency, and t c \ can be approxi- 
mated by the mean time between collisions described by 
the classical Boltzmann collision integral 622- In static 
equilibrium, this collision rate for a uniform gas is given 
by 
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(4) 



where n is the density of atoms, a = 8ira 2 is the quan- 
tum collision cross section, and v is the average speed of 
an atom in the gas. Both above and below TbeCj the 
analogous collision time corresponding to collision pro- 
cesses described by C22 in (H) will give an estimate of 
the lifetime of a single-particle excitation in the ther- 
mal cloud. This is distinct from the physics given by 
C12, which describes the collisions of condensate atoms 
with atoms from the thermal cloud. In particular, one 



finds that if the thermal cloud is described by the equi- 
librium Bose distribution (/ = /°), then C 2 2[/°, c f ) ] = 
but Ci2[/°,$] 7^ 0- Thus C12 will give rise to damping 
of condensate oscillations even when the thermal cloud is 
treated statically. Of course, in the collisionless region, 
there is another source of damping arising from the dy- 
namical mean-field coupling between the condensate and 
thermal cloud that is also included in (|l|) and (||) ; this is 
Landau damping , which will be discussed below. 



A. Static Popov approximation 

In the present work, we use these equations to calculate 
the damped normal modes of the condensate given by 
the solutions of ([]]) assuming that the non-condensate 
atoms always remain in static thermal equilibrium. For 
our model, this means we take 



/(p,r,i)~/°(p,r) 
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e (3[p 2 /2m+U (r)-p,o] _ 1 



(5) 



where fio is the equilibrium chemical potential of the non- 
condensate and Uo(r) = U cxt (r) + 2g[n c o(r)+ho(r)]. The 
detailed analysis given by ZNG shows that the Bose- 
Einstein distribution in (|^) is a stationary solution to 
(Q) when the condensate and non-condensate are in dif- 
fusive equilibrium, which requires jlo — /i c o, where /i c o is 
the equilibrium chemical potential of the condensate as 
described by ([j]). 

Using our finite T "static Popov" approximation, equa- 
tion (fy) can be simplified to 
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which describes the condensate motion within the static 
thermal cloud. Here no is the equilibrium density of the 
non-condensate and the damping term Rq is calculated 
using 

r? 2 M) = J j0^C 12 [f°(p,r)Mr,t)]. (7) 

Notice that r° 2 (r,t) depends on time now only through 
$(r,t). Using the explicit general expression for C12 
given in Eq. (23b) of ZGN', one finds (see also Ref. j§) 

r?2(r ' t) = ^| [e ^ ( ' i0 " C(r ' t)) " l] ' (8) 
where we have defined the C12 collision time 

— 7~~l\ = ^ 2g .^ + 7 [ ^Pi / d Pz Ap3<*(Pc + Pi - P2 - Pa) 
n 2 (r,t) (2tt) 5 ?)/ J J J 

x 6(e c + e Pl ~ e P2 - e P3 )(l + ft) ft ft . (9) 

Here the condensate atom local energy is e c (r, t) = 
/i c (r, t) + \mv 1 (r, t) with the non-equilibrium condensate 
chemical potential 
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Mc(r,t) = - 



2m^/n 



(10) 



The condensate atom momentum is p c = tov c , and 
fi = /°( r j Pi)- We have introduced the usual condensate 
velocity defined in terms of the phase of the condensate 
$(r,t) = y/n c {r, t)e ie ^ as v c = hV0(r,t)/m. A closed 
set of equations for $(r, t) is given by @ and its complex 
conjugate combined with (||) - (|Io|). 

We note that in terms of n c and v c , (Q) is completely 
equivalent to the coupled equations [|| 



^ + v-(n c v c ) = -r; 2 [/°,$] 



(9f 



V v c = - V^ c . 



(11) 



It is easy to see from (||) that when the condensate 
is in equilibrium with the thermal cloud according to 
Me — ► Mco = Aoi ri2( r ;^) then vanishes. It is clear that 
the description of the system given by equations (||) and 
(Q), or equivalently (pi]), is valid only if the condensate is 
slightly perturbed from equilibrium and the condensate 
motion is essentially uncoupled from that of the ther- 
mal cloud, which we can then treat statically. In or- 
der to describe the condensate oscillations about equi- 
librium, we use the quantum hydrodynamic variables 
n c (r, t) = n c0 (r)+Sn c (r, t) and v c (r, t) = 5v c (r, t), where 
n c0 (r) is the equilibrium density of the condensate with 
the associated equilibrium chemical potential fi C Q. Alter- 
natively, one may work with the fluctuations of <I>(r,i) 
and derive coupled Bogoliubov equations |2l],|23|,Q gen- 
eralized to include the effect of the Rq damping term. 
This generalization will be discussed elsewhere |22] . 



B. Finite-T Stringari wave equation 



From (11), we can obtain linearized equations of mo- 
tion for the condensate fluctuations Sn c and <5v c . We use 
the fact that, to lowest order in the fluctuations from 
static equilibrium, (ph reduces to 
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T 12\ r ) 
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where the "equilibrium" C12 collision rate is defined by 

1 2ff 2 
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gn^il + f^fVl (13) 



In the present paper, we restrict ourselves to the Thomas- 
Fermi limit, valid for large N c , 
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V • (n c0 <5v c 
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Sn c 



m— — = -gVdn c . 
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(14) 
(15) 



The collision time r'(r) describes collisions between the 
condensate and non-condensate atoms when the conden- 
sate is perturbed away from equilibrium, 
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3"co(r) 1 



r'(r) k B T 7f 2 (r) 



(16) 



In the TF approximation the equilibrium distribution re- 
duces to /P = [ e /3(pf/2™+snco) _ The new term on 
the right-hand side of ( |l4| ) causes damping of the con- 
densate fluctuations due to the lack of collisional detailed 
balance between the condensate and the static thermal 
cloud. We note that this collision time is only a function 
of r through its dependence on the static condensate den- 
sity n c0 (r). Plots of l/r'(r) will be discussed below. 

We can easily combine (|l4|) and ( |l5| ) to obtain what 
we shall refer to as the finite T Stringari wave equation 



d 2 8n c 
dt 2 
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V • (n c0 VSn c ) 



1 dSn c 



(17) 



Equation ( |l7j ) is the main result of this paper. If we ne- 
glect the right-hand side, we obtain the undamped finite 
T Stringari normal modes 5n c (r,t) — 5ni(r)e~ tuJit given 
by the solution of || 



-^V • [n c o(r)V<Srci(r)] = uf5m(r). 
m 



(18) 



As has been noted by several authors in recent 
papers Jll| , ^2| . p3| , n c0 (r) at finite T can be well approx- 
imated by the TF condensate profile at T = but with 
the number of atoms in the condensate N C (T) now being 
a function of temperature, since the static mean field of 
the non-condensate plays such a minor role. With this 
approximation for n o(r), the solutions of the finite T 
Stringari equation ( |l8| ) will be identical to those at T = 0, 
since the T = Stringari frequencies do not depend on 
the magnitude of A^ c . Of course, as shown by calculations 
solving the coupled Bogoliubov equations K,^3| , the TF 
approximation breaks down when N c < Kr. Thus the 
condensate collective mode frequencies will always be- 
come temperature dependent close to TbeCi where the 
TF approximation is no longer valid. We generalize our 
present discussion to deal with this region in Ref . @] . 

We can use the undamped Stringari modes as a basis 
set to solve ( |l7| ) and find the damping of these modes. 
Writing 5n c (r) = J2i c i^ n i( r )^ an d using the orthogo- 
nality condition j drSni(r)5rij(r) 
following algebraic equations for the coefficients c 



Sij, one obtains the 



LO Ci — LO i Ci~ ILU 



where 



jij = / drSni(r)5nj(r)/T'(r) . 



(19) 



(20) 
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Assuming the damping is small (we are in the collisionless 
region) , ([l9]) is easily solved using perturbation theory by 
setting 7y = for i ^ j . This gives the damped Stringari 
frequency (to lowest order) Qi = ioi — iTi, with 



2*1 = - [ dr- Ui ^ 



r'(r) 



(21) 



This result for Ti is reasonable, namely it involves an 
average over l/r'(r) weighted with respect to the un- 
damped density fluctuations of the Stringari wave equa- 
tion ( p!s|) . We find that coupling to other modes (7^ 7^ 0) 
is extremely small. 



III. RESULTS 
A. Homogeneous gas 

Before treating the trapped gas, it is useful to first ap- 
ply our theory to a homogeneous gas, which was con- 
sidered previously in Ref. p5fl in connection with the 
collision-dominated hydrodynamic region. For a homo- 
geneous gas, t' is independent of position and then j2l| ) 
reduces to I\ = 1 /2r'. Although our model in the present 
paper applies only to the collisionless region, it is useful 
to compare the inter-component collision time in both 
the collisionless and hydrodynamic regimes. In ZNG, it 
was shown that the inter-component collision time 7), in 
the hydrodynamic region is given by 7), = err', where the 
temperature-dependent factor a (not to be confused with 
the collision cross section) depends on various thermody- 
namic functions. In Fig. 1 we compare 1/r' and 1/t m 
as functions of T. We see that a dramatically alters the 
inter-component relaxation rate l/r M appropriate to the 
hydrodynamic regime, as compared to 1/r' involved in 
the collisionless regime. For completeness, in Fig. 1 we 
also plot the often-used classical collision time given by 
(@) as well as t® 2 defined in (^|). 

In a uniform Bose gas at finite temperatures, the Lan- 
dau damping (u> — cq — iT^) of condensate modes has 
been evaluated in several recent papers [|s|-|l4| . Working 
within the full second-order Beliaev approximation, one 
finds 



3it\ ak^,Tq 
8~ 



(22) 



This is clearly quite different from our inter-component 
damping V = l/2r', as plotted in Fig. 1. Landau damp- 
ing originates from the interaction of a condensate collec- 
tive mode with the excitations of the thermal cloud but 
is not associated with C12 collisions which give rise to t'. 




T/T 

c 

FIG. 1. Collision rates in a homogeneous Bose-condensed 
gas, normalized to the collision rate of a classical gas at 
the BEC transition temperature t c Y 1 (Tbec)- We have taken 



gn — O.lfcBTBEC, where n is the total density. 
Refs. ElEi. 



See also 



In the context of our generalized GP equation in (Q), 
Landau damping comes from the fluctuations in the ther- 
mal cloud induced by the condensate mean field, 



5n = x (2gSn c 



(23) 



In the finite temperature region of interest, xo can be 
approximated as the density response function of a non- 
interacting gas of atoms with a spectrum e p and chemical 
potential [i c q. For a uniform gas, one sees that using j23| ) 
m (0), with R = 0, gives condensate modes satisfying 



c q (1 + igxa) and thus lo = cq — iT^, where 
T L = 2gcqlmxo(q, uj = cq) . 



(24) 



Evaluating Imxo in the limit of small q Q, one finds 
Tl = ^ak^Tq/U. Apart from the slightly larger numer- 
ical coefficient, this agrees with the exact result given in 
equation (||) ||o). 



B. Trapped gas 

We now turn to explicit calculations of the inter- 
component damping rate using our model for a trapped 
gas. In order to calculate Ti for a given mode, the 
equilibrium chemical potential fi c o = /to must be cal- 
culated self-consistently for a given total number N of 
atoms at a given temperature T. In the TF approxima- 
tion, the procedure is straight forward (see, for exam- 
ple, Ref. ^7j). In the following, we consider a harmonic 
trap with axial symmetry C/ oxt (r) = ^mu>p(p 2 + X 2 z 2 ), 
where A = lu z /lu p is the anisotropy parameter. In the 
TF approximation, the condensate density takes the ex- 
plicit form n C Q = [fi C Q — muip(p 2 + X 2 z 2 )/2]/g within 
the TF radius, and the condensate chemical potential is 
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Mco = 5^p[15AiV c a//9o] 2 / 5 , where po — \fhJmuT p . The 
form of the Stringari normal modes Srii(r) is given ex- 
plicitly in the literature H^lJ- We mainly consider the 
breathing mode [n = 1 ,1 = 0) for which wio = V5uj p , for 
A = 1. 

We choose experimentally accessible parameters in the 
following calculations for the collisionlcss region. How- 
ever, we do not compare our results to the two available 
experiments on damping of normal modes at finite T, 
since the TF approximation is not valid for most of the 
JIL A data (2£| , and the MIT experiment |^] is approach- 
ing the collision-dominated hydrodynamic limit where 
the dynamics of the condensate and non-condensate be- 
come more strongly coupled. For 87 Rb the scattering 
length is a ~ 5.7 nm We first consider a spherically 
symmetric trap A = 1, with trap frequency v r = 10 Hz, 
and we take N = 2 x 10 6 . In the collisionlcss limit, we 
require 0JiT c \ ^> 1, taking t c \ as defined in ([|). For a 
trapped gas, we obtain an upper limit on 1/t c \ by tak- 
ing the density in the center of the trap n(r = 0), which 
gives 1/t c i = %a 2 N ui^m / {irk-sT) . For the parameters we 
use, wiot c i ss 19 (compared to wo2T"ci w 20 for the JILA 
data (28), and w 02 t c1 w 2 for the MIT data @). 



havior of l/r'(r) just seems to be mimicking the behavior 
of the non-condensate density n(r), which we plot in the 
inset of Fig. 2b along with the condensate density. The 
condensate mean-field pushes the non-condensate density 
out of the center of the trap, a well known result pi) . 
We also show the breathing mode density fluctuation in 
Fig. 2b. The sharp cusp of n(r) and the sudden drop 
of 5rii{v) and l/r'(r) at the condensate boundary are 
all unphysical artifacts of the TF approximation. In- 
clusion of the kinetic energy pressure in a more accu- 
rate calculation would have the effect of smoothing out 
this behavior at the boundary. To illustrate the effect of 
the kinetic energy pressure, we also show in Fig. 2b the 
breathing mode obtained by solving the T = coupled 
Bogoliubov equations. We estimate that an improved 
treatment which includes the kinetic energy pressure will 
modify our estimate of Tio by about 10 — 20% 

In Fig. 3a we plot the damping rate Tio for the breath- 
ing mode (n = 1, I = 0) shown in Fig. 2b as a function 
of temperature up to T = 0.95T B ec, where jV c «7x 10 4 . 
At higher temperatures, the Thomas-Fermi approxima- 
tion will start to break down and the mode frequencies 



become temperature dependent [24 2 




0.06 



r/R Tp 

FIG. 2. Positional dependence of various quantities. In (a) 
we plot l/r'(r) normalized by its value at the TF radius -Rtf- 
In (b) we show the density fluctuation Smo of the Stringari 
breathing mode (solid) for a spherically symmetric trap. We 
also show the exact T — Bogoliubov mode (dashed) for 
N C (T = 0.9T B ec) = 2.3 x 10 5 . Both solutions are normalized 
to unity, J 8n\ (v)dv — 1. The densities of the condensate 
and thermal cloud are plotted in the inset for T = 0.9Tbec- 

In Fig. 2a we plot l/r'(r) vs. position for T = 0.9Tbec 
and T — 0.5Tbec- We see that the collision rate in- 
creases steadily up to the condensate boundary, but as T 
increases, l/r'(r) becomes relatively constant. The be- 
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FIG. 3. Normal mode damping rates vs. temperature. 
In these plots, the damping rates are normalized by their 
corresponding mode frequencies and we only plot up to 
T = 0.95Tbec, above which the Thomas-Fermi approxima- 
tion will start to break down. In (a) we show the damping rate 
for the breathing mode (n = 1, I = 0) of a spherically symmet- 
ric trap, where the solid line corresponds to inter-component 
collisional damping given in (pl|). In (b) we show damping 
rates for the quadrupole modes (n = 0, I = 2) in a cylindrical 
trap. The solid line is for the m — mode and the dot-dashed 
line is for m — 2. 

Landau damping of condensate modes in trapped gases 
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has also been discussed in some detail in recent pa- 
pers 

@@- These 

papers give results that are in quali- 
tative agreement with the expression for a uniform gas in 
(p2|), with q — ujp/c and evaluating the Bogoliubov sound 
velocity c for the density at the center of the trap ||. 
This simple estimate for Landau damping is plotted in 
Fig. 3 for comparison. We note that it is larger but com- 
parable to the inter-component collisional damping that 
we consider. Clearly, a fully satisfactory theory of finite 
T damping of normal modes must include both Landau 
damping as well as the damping we consider in this paper 
due to the condensate being out of diffusive equilibrium 
with the non-condensate. 

Our theory is easily applied to anisotropic traps. In 
Fig. 3b we show the damping of m = 0, 2 quadrupole 
modes for an axially symmetric trap with A = y/8. Here 
we choose a slightly tighter trap v r = 23 Hz, and we take 
N = 1 x 10 6 (in this case, N c « 3 x 10 4 at T = 0.95T B ec)- 
For these parameters, we find u)2ot c \ ~ 6. In Fig. 3b we 
see that Landau damping is about twice as large as our 
inter-component collisional damping. 




10 6 N 

FIG. 4. Damping rate vs. T and total TV of the breath- 
ing mode. Here we plot the surface of I\o and its con- 
tours projected onto the plane below. The solid white line at 
N — 2 x 10 6 corresponds to the solid line plotted in Fig. 3a. 
The longer line illustrates that in experiments, TV decreases as 
the temperature is lowered due to evaporative cooling. The 
upper edge of the surface corresponds to the critical temper- 
ature (kBTBEc/tiujp) = 0.947V 1 / 3 , above which Vio vanishes. 

It is instructive to also consider the dependence of the 
mode damping I\ on the total population N. In Fig. 4, 
we show a shaded surface plot of Tio for the breathing 
mode in an isotropic trap as a function of T and N. 
The white line at N — 2 x 10 6 corresponds to the solid 
line plotted in Fig. 3a. As one might expect, the inter- 
component damping rate increases with increasing total 
population (since the density is increasing). It is also 



important to realize that in current experiments, the data 
taken is for a broad range of N due to evaporative cooling 
losses p^ , p9| . The idealized fixed- TV line can never be 
achieved in practice and one is instead dealing with a 
curve like the longer line on the surface in Fig. 4. 

IV. CONCLUSION 

In summary, we have calculated a new damping mech- 
anism of condensate collective modes due to collisions 
with the thermal cloud, based on the finite-T equations 
derived in Ref. 0. The essential mechanism involves 
the lack of diffusive equilibrium between the conden- 
sate and the thermal cloud Q, which also plays a key 
role in the theory of condensate growth fl9| , po|| . Here 
we have carried-out the first explicit calculation of this 
damping mechanism for a trapped gas in the collision- 
less regime (this inter-component damping has recently 
also been evaluated in the collision-dominated hydrody- 
namic regime p2[). In recent discussions of the damping 
of condensate collective modes in the collisionless region, 
the mechanism we consider is omitted. One instead fo- 
cuses on the dynamical mean-field coupling between the 
condensate and thermal cloud, which gives rise to Landau 
and Beliaev damping 0-|l4|. While we have not consid- 
ered it in detail, we have indicated how we could include 
Landau damping by considering the non-condensate fluc- 
tuations in (Q) induced by the condensate mean field . 
Comparing Landau damping to the additional mecha- 
nism we have calculated, we find that the two are com- 
parable in size. Further experimental studies of the col- 
lective modes at finite temperatures are needed to clar- 
ify the relative importance of these different sources of 
damping. 

In this paper, we have argued that a good first estimate 
of the inter-component damping of condensate collective 
modes can be obtained by coupling it to a static thermal 
cloud. A more systematic theory is clearly desirable in 
which the collisionless dynamics of the thermal cloud is 
allowed for. However, as noted in the introduction, we 
do not believe that this will lead to significant corrections 
to the inter-component damping of out-of-phase conden- 
sate modes in which the motion of the thermal cloud is 
not significant. In future work, we hope to discuss the 
damping due to Cn collisions of collective modes that 
mainly involve the motion of the thermal cloud, with the 
condensate being treated statically. 

This work grew out of discussions with E. Zaremba 
about the role of Cn collisions in the collisionless region. 
We also thank T. Nikuni for useful comments. This re- 
search was supported by a grant from NSERC. 
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